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Using a classical equation of motion, dynamics of the phase is analyzed in the Inductive 
Single-Electron Transistor (L-SET) which is a promising new system suitable for quantum 
measurement with ultimate sensitivity and low back-action. In a regime of nonlinear dynamics, 
a shift of the oscillator resonant frequency is discovered which has a direct analogy to the 
switching of a dc-biased Josephson junction into voltage state. Results are reviewed for the 
predicted charge sensitivity, and it is shown that a performance challenging the best rf-SETs 
is foreseeable with the new device. 

1 Introduction 

The Radio-Frequency Single-Electron Transistor (rf-SET) has proven a suitable workhorse for 
high-frequency charge measurements. Basically, no other device has been able to track dynamic 
single-charge transport at MHz frequencies, a property relevant especially from the point of 
view of characterization and readout of superconducting qubitsl^. However, there are reasons to 
search for an alternative technology to overcome limitations principally due to the dissipative 
nature of the rf-SET"^; (1) it is not a truly quantum-limited detector, (2) relatively high power 
dissipation required for operation heat up the sample and surrounding qubits considerably, (3) 
it would be beneficial to be able to combine qubit and detector into a single device. 

Purpose of the present paper is to present an intuitive picture of the operation a new kind of 
radio-frequency charge detector, the Inductive Single-Electron Transistor (L-SET), introduced 
recently by Sillanpaa et al. '^ both theoretically and experimentally. The L-SET is based on 
reactive readout of the Josephson inductance of a superconducting SET (SSET) in a resonator 
configuration. The schematics is ultimately a non-dissipative, high-bandwidth, and sensitive 
electrometer with the important property that it lacks the shot noise and thus excessive back- 
action inherent in dissipative detectors due to sequential tunneling. 

2 The L-SET circuit 

We use the circuit shown in Fig. ^ where the SSET is coupled in parallel to an LC oscillator 
resonant at the frequency /o = l/(27r)(LC)~^'^. The resonance of the whole system, fp = 
l/(27r)(L II LjC)^^/"^ > /o depends on the SSET effective Josephson inductance {Lj)~^ = 
(27r/$o)^-^j in parallel with the external inductor L. Here, $0 = ^/(2e) is the flux quantum. The 
effective values of Josephson coupling Ej = d'^E{qg,(l))/d(j)'^ and critical current Iq ~ 2eE*j/% 
(not exact due non-sinusoidal current-phase relation) are due to the quantum band structure 
E{qg,(p), where (j) is the phase difference (p = 2e/h ^^ V{t)(it across the SSET. Since L*j{qg) can 
have a substantial dependence on the gate charge qg = CgVg, the value of the resonant frequency 
fp can be used for building a sensitive detector. Its bandwidth A/ ~ fo/Qe, where Qe is the 
coupled quality factor, is typically in the range of tens of MHz. 

3 Phase dynamics 

Let us consider the SSET as a single Josephson junction which has the current-phase relation 
I J ~ /q sin((;/)). Using the RCSJ-model, a classical equation of motion can be derived for the 
phase difference (p across the resonator. The capacitance is due to the shunting C (Fig. P) and by 
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Figure 1: Schematics of the L-SET electronics. The component values refer to sample 2. The resonance is read 
by measuring amplitude and phase of the reflection coefBcient T = {Z ~ Zo)/{Z + Zq), where Zq — 50fl, and 

Z{L*j) the resonator impedance (including Cc)- 



the resistor R we model the small residual dissipation. The resulting dynamics has an intuitive 
mechanical analog. 

Current / in Fig. ^ is divided in the junction, R and C, 

/ = /q sin(0) + V/R + CV = /q sin((/)) + h/{2eR)<j) + Ch/{2e)'^. 
The system is driven trough the coupling capacitor Cc- The originating current lext is divided, 
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where Idrc is the current circulating in the loop through the external inductor 
L. It is simply related to (p via the magnetic flux in L, <1> = ^o/2tt(J) = Lldrc- 

Combining, we have classical non-linear equation of motion for internal dynamics of the 
phase 
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which is analogous to a mechanical particle moving in a sinusoidally modulated parabolic potential 
V{(p) = ^Q/{87r^L)(p'^ — Ej cos(0) and subject to a force lext- Let us take a drive lext = le cos{ujt), 
and assume that the response is dominated by components at the drive angular frequency, 
(p = xsin(a;t) + ycos{ujt). By keeping terms of the zeroth and the first harmonic, non-analytic 
equations are derived for the quadrature amplitudes x and y, 



-Mij^x - cuoy + kx + 21^ [Ji{x)Jo{y) + Ji{x)J2{y)] = 
-Mio'^y + CLOX + ky-Ie + 21^ [Jo{x)Jiiy) + Ji(y) J2(x)] = 



(2) 
(3) 



where Jn are Bessel functions of the first kind and order n, M = Ch/{2e), c = h/{2eR), and 
k = $o/(2vrL). In the absence of friction, R -^ oo, the response is in phase with the drive, 
X = 0, and Eqs. 2, |31 simplify to 



MuJ^y + ky-Ie + 2I*Ji{y)=0 



(4) 



which allows for a simple graphical interpretation. 

The solution for the amplitude y is the intersection of the line z = {k — Muj'^)y — Ig, and 
the Bessel term z = — 2/9 Ji(2/). Slope of the line, k — Muj"^, decreases with increasing drive 
frequency oj. At the resonance frequency ujji{Ie), the phase of the response shifts by 2tt, which 
in this picture corresponds to change of sign of the solution y. While increasing the drive 
frequency (see Fig. |2I), the solid lines a..-d, plotted for increasing frequency but constant drive 
le = 0.5/q, rotate clockwise about the point z = —le- At relatively low drive, the favored 
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Figure 2: Eq. 2] illustrated, with M = k = Iq — 1, and o^o = ^Jk/M = 1. The smallest-amplitude solution is 
marked by hollow symbols, and the resonance condition by solid symbols. For lines e and g, ui = 1.0 and 1.05, 

respectively. 
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Table 1: Parameters of the samples and their tank LC oscillators discussed in the text. Ej — hA/{8e'^^RsET) is 
the single-junction Josephson energy, and Ec = e^ /{2Cj:) is the charging energy. 

solution which has the smallest absolute value and hence the lowest energy, changes sign at 
the frequency lor ~ 1.19a;o where the line d tangents the Bessel function at the point marked 
with solid circle. We note that the line d has a clearly negative slope, which corresponds to a 
relatively large value of ujr. 

While increasing the drive, at the particular value le — 1-2/q the sign of the smallest solution 
changes when the line is horizontal, as depicted by the line e, which corresponds to a smaller 
frequency than in the previous case. At this to, the solution changes from y ~ 1.8 (solid diamond) 
to y = — cxD, and approaches again zero (line g) at higher frequency. The infinity is due to absence 
of dissipation. At still higher drive, the resonance frequency does not change in this simplified 
picture. Note also that if the Bessel term due to the Josephson contribution -^ 0, the resonance 
condition corresponds to a horizontal line as well. 

We re- formulate the statements: At the "critical drive" le — l-2/o the resonance frequency 
cvji shifts from a high cj/j^ = 27r/p which includes the Josephson-contribution, to the lower 
ujji = 27r/o which is determined only by L and C. 

This shift of resonance frequency of the L-SET oscillator at a certain "critical power" Pc of 
the microwave drive is in a direct analogy to the switching of a DC-biased Josephson junction 
into voltage state ". Of course, a similar phenomenon should take place equally well in a single 
junction than the SSET. 



4 Experimental results 

Results of two samples (see Table P) are discussed. We focus here on the phase dynamics of 
sample 1. Since its Ec was comparable to temperature, it did not work as a good electrometer, 
however, it was therefore rather close to a classical junction supposed in the calculations. Sample 
2 (discussed in Ref.™, see also section|SJ was a sensitive detector, with a gate shift of the resonance 
frequency of 15 MHz. 

In Fig. 121 (a) we plot the frequency response for sample 1 as a contour, a suitable repre- 
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Figure 3: Frequency response (|r|) of sample 1 as a function of increasing drive amplitude, (a) experiment; (b) 
numerical calculation using Eqs. 2,|3]and circuit relations for Z. The gray scales in (a) and (b) are scaled from 1 

to the smallest |r| in each case. 

sentation for the dynamics as a function of drive strength. The shift of resonance frequency 
from the plasma resonance at /p = 710 MHz to the tank resonance at /o = 610 MHz is clearly 
visible at a critical power Pc — —102 dBm. A prediction from our formalism which uses only 
independently determined parameter values, Fig. 13] (b), agrees qualitatively with the features of 
the experimental data (yet the resonance depth is different due to probably too small value of 



5 Charge sensitivity 

The bottleneck for detector sensitivity in the present schematics is the preamplifier having a noise 
temperature of 5 K. For sample 2 in the discussed operation mode of small plasma oscillations, 
we measured a charge sensitivity Sg ~ 2 x 10~^e/vHz. At higher drive which corresponds to 
several periods of the cos(0) Josephson term in the potential ("anharmonic mode"), such that the 
transduction is not due to tuning of /p but rather, due to tuning of non- linear phase dynamics, 
a better sensitivity of 1.4 x 10~^e/vHz was achieved. Since the sensitivity is predicted to 
improve roughly as {Ej/Ec)~^, the record numbers of the rf-SET, Sg ~ 3 x 10~^e/\/Hz, should 
be beatable at experimentally accessible values of Ej/Ec ~ 0.4 in the anharmonic mode, or 
with Ej/Ec ~ 0.15 in the plasma oscillation mode, both using an Al SSET in the present 
experimental scheme'^. 
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